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KINEMATICS OF THE X-SHAPED MILKY WAY BULGE: 

EXPECTATIONS FROM A SELF-CONSISTENT A-BODY MODEL 

Yujing Juntai Shen^ ^ Zhao-Yu lP, Shude Mao^’"^’^, Martin C. SmithP 

R. Michael Rich Andrea Kunder^, and Chao Liu^ 

ABSTRACT 

We explore the kinematics (both the radial velocity and the proper motion) of the vertical X-shaped feature 
in the Milky Way with an A-body bar/bulge model. From the solar perspective, the distance distribution of 
particles is double-peaked in fields passing through the X-shape. The separation and amplitude ratio between 
the two peaks qualitatively match the observed trends towards the Galactic bulge. We confirm clear signatures 
of cylindrical rotation in the pattern of mean radial velocity across the bar/bulge region. We also find possible 
imprints of coherent orbital motion inside the bar structure in the radial velocity distribution along 1 = 0°, where 
the near and far sides of the bar/bulge show excesses of approaching and receding particles. The coherent orbital 
motion is also reflected in the slight displacement of the zero-velocity-line in the mean radial velocity, and the 
displacement of the maximum/minimum in the mean longitudinal proper motion across the bulge region. We 
find some degree of anisotropy in the stellar velocity within the X-shape, but the underlying orbital family of 
the X-shape cannot be clearly distinguished. Two potential applications of the X-shape in previous literature 
are tested, i.e., bulge rotation and Galactic center measurements. We find that the proper motion difference 
between the two sides of the X-shape can be used to estimate the mean azimuthal streaming motion of the 
bulge, but not the pattern speed of the bar. We also demonstrate that the Galactic center can be located with the 
X-shape, but the accuracy depends on the fitting scheme, the number of fields, and their latitudinal coverage. 
Subject headings: Galaxy: bulge - Galaxy: kinematics and dynamics - Galaxy: structure 


1. INTRODUCTION 

The Galactic bulge holds important information towards 
understanding the formation and evolution of the Milky Way 
(MW). However, our edge-on perspective in the disk plane 
and the large and variable dust extinction especially in op¬ 
tical bands make it a non-trivial task to study the MW bulge. 
Despite these difficulties, our knowledge of the Galactic bulge 
has advanced greatly in the past two decades. Numerous stud¬ 
ies indicate that the MW bulge contains a bar structure. For 
example, the non-circular motions of gas in the inner Galaxy 
was explained with bar perturbation (de Vaucouleurs 1964; 
Fux 1999; Englmaier & Gerhard 1999; Bissantz & Gerhard 
2002). Likewise, the boxy shape of the MW bulge revealed 
by COBE/DIRBE photometry (Dwek et al. 1995) was charac¬ 
terized as the perspective effect of a tilted bar (Blitz & Spergel 
1991; Binney et al. 1997). More evidence for the barred mor¬ 
phology emerges from star counts using distance indicators 
like Red Clump (RC) stars (Stanek et al. 1994; Rattenbury 
et al. 2007; Cao et al. 2013). Kinematic information of 
the MW bulge stellar populations provides additional insights 
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into its origin. The Bulge Radial Velocity Assay (BRAVA) 
revealed the fast cylindrical rotation of the Galactic bulge ex¬ 
tending up to IZ?I = 8°, which is inconsistent with the classi¬ 
cal scenario that the bulge is built in mergers and dominated 
by random motions (Howard et al. 2008, 2009; Kunder et 
al. 2012). Shen et al. (2010) demonstrated that a pure-disk 
galaxy with a spontaneously formed bar structure can match 
the BRAVA observation of the radial velocity and velocity dis¬ 
persion towards the bulge region well. They also found that 
the mass of a classical bulge, if in existence, may be relatively 
small. Therefore, the Galactic bulge is likely built through 
bar-induced internal dynamical processes, with minor contri¬ 
bution from a merger-built classical bulge. 

Recently, some studies suggested that the Galactic bulge 
shows the signature of a vertical X-shape. Towards the bulge 
region, the apparent magnitudes of RCs in certain fields show 
double-peaked distributions, with increasing peak separations 
towards higher latitudes (McWilliam & Zoccali 2010; Nataf 
et al. 2010; Saito et al. 2011; Ness et al. 2012; Nataf et 
al. 2015). Since RCs are well-calibrated standard candles 
with a roughly constant absolute magnitude, such a double- 
peaked feature seems to indicate the existence of a vertical 
X-shaped feature in the Galactic bulge. More detailed stud¬ 
ies revealed that certain latitudinal slices of the Galactic bulge 
have two density lumps aligned with the bar, which indeed re¬ 
sults in the observed double-peaked apparent magnitude dis¬ 
tributions of RCs (Saito et al. 2011). Based on the RCs in the 
VISTA Variables in the Via Lactea (VVV) Survey, the sym¬ 
metrized three-dimensional bulge model of Wegg & Gerhard 
(2013) also clearly illustrates the X-shape in its side-on view 
and surface density slices. Further Made-to-Measure mod¬ 
elling reveals an off-centered X-shape in the Galactic bulge, 
and the peanut/X-shaped component accounts for over 20% 
of the bulge stellar mass (Portail et al. 2015a). The X-shaped 
bulge morphology is frequently seen in edge-on disk galaxies 



2 


(e.g. Whitemore & Bell 1988; Lutticke et al. 2000; Bureau 
et al. 2006), and is presumably similar to the structure in our 
own Galaxy. 

Using a simple and self-consistent A^-body simulation, Li 
& Shen (2012) showed that such an X-shape can arise nat¬ 
urally during the evolution of the Galactic bar. Numerical 
simulations demonstrated that a bar forms via the bar instabil¬ 
ity, and then thickens vertically due to the buckling instability 
(Combes & Sanders 1981; Raha et al. 1991). The vertical 
motion in a peanut/X-shape bulge is enhanced after the buck¬ 
ling episode. Some studies proposed that xl orbits with 2:1 
vertical Lindblad resonance, namely banana orbits for their 
side-on appearance, can contribute to the peanut/X-shaped 
feature in edge-on galaxies (Pfenniger & Friedli 1991; Pat- 
sis et al. 2002; Athanassoula 2005). Great efforts have been 
made to understand the roles of banana orbits (and their as¬ 
sociated orbits) in forming the peanut/X-shaped bulges (e.g. 
Patsis & Katsanikas 2014a, b). However, recent work by Por- 
tail et al. (2015b) classified orbital families in the peanut/X- 
shaped bulges, and suggested the brezel-like orbit as the main 
contributor to the X-shape. Since all of these orbital fami¬ 
lies can potentially support the peanut/X-shaped feature, ad¬ 
ditional kinematic information is required to distinguish their 
contributions in the Galactic X-shape. If the Galactic X-shape 
is dominated by certain orbital families, its stellar kinematics 
is expected to refiect the underlying orbital motion. Thus, the 
stellar kinematics is essential to better understand the Galactic 
X-shape. 

Currently, only a handful of observational studies have ex¬ 
plored the kinematics inside the X-shape. Vasquez et al. 
(2013) (hereafter VI3) analyzed the stellar motion in the 
X-shape using RCs in (/,/?) = (0°,-6°). With the three- 
dimensional velocity available for 352 stars, they found weak 
anti-correlation between the longitudinal proper motion and 
the radial velocity in both bright and faint RCs. Meanwhile, 
no significant correlation was seen between the latitudinal 
proper motion and the radial velocity in their sample. They 
interpreted the correlation with stars on elongated bar orbits, 
possibly the streaming motion along the arms of the X-shape. 
Poleski et al. (2013) (hereafter PI3) investigated the proper 
motion of RCs in fields at Z? ^ 5°. Between the probability- 
distinguished closer and further arms of the X-shape, the 
mean proper motion difference is asymmetric in both I and 
b directions, which is linear for -0.1° < / < 0.5°, but roughly 
constant for -0.8° < / < -0.1°. They attributed the linear side 
to the streaming motion within the X-shape. From a numer¬ 
ical perspective, Gardner et al. (2014) (hereafter G14) first 
explored the kinematic properties of X-shaped bulges. They 
compared the kinematics at the closer and further sides of the 
bar/bulge, and suggested that X-shaped bugles show coherent 
signature of a minimum along / = 0° in the mean radial veloc¬ 
ity difference between the closer and further sides. Although 
efforts have been made in both observations and simulations, 
the kinematics of the X-shape is still poorly understood. 

In this paper, we exploit the power of simulations to present 
a closer look at the X-shape in our own Galaxy. Using an N- 
body model of the Milky Way-like barred galaxy, we explore 
the kinematic features of the Galactic X-shape, and compare 
our model to the observations and other simulation works. 
The numerical model here can not represent the actual Galac¬ 
tic bulge in all aspects, but its self-consistency and simplicity 
still reveal the general properties of the X-shaped bulge in 
such exploration and comparison. Our results lead to a bet¬ 
ter understanding of the existing observations, serve as an in¬ 


dependent check of other numerical works, and make predic¬ 
tions for future observations. In Section 2, we briefiy describe 
the A-body model used in this paper. Results and comparisons 
with existing observations are shown in Section 3. In Section 
4, we discuss the nature of the X-shape by analyzing the cor¬ 
relations between proper motions and radial velocities, and 
test two potential applications to derive Galactic parameters 
based on the Galactic X-shape. Our results are summarized in 
Section 5. 

2. MODEL 

In this work, we use the same A-body model as in Shen 
et al. (2010) and Li & Shen (2012). Initially it contains a 
cold (Toomre’s 1.2) exponential stellar disk of 10^ par¬ 
ticles rotating in a rigid dark matter halo potential. The bar 
forms rapidly within the first gigayear, and becomes thick¬ 
ened vertically in the inner region via the buckling instability. 
After ^2.4 Gyr, the structure of this model reaches a quasi¬ 
steady state, with a bar pattern speed of ^ 40 km s“^ kpc“^ 
(Shen 2014; Molloy et al. 2015). The snapshot at ^ 4.8 Gyr 
was analyzed in Shen et al. (2010) and Li & Shen (2012). 

Assuming a distance of the Sun to the Galactic center 
(GC) (Ro ^ 8.5 kpc) and a 20° angle of Sun-GC line to the 
bar major axis, this model yields excellent matches with the 
mean radial velocity and velocity dispersion measurements 
of BRAVA. The X-shape in this model was highlighted after 
subtracting a smooth component from the edge-on projection 
(Li & Shen 2012). The distance distributions in the model 
towards the bulge fields can qualitatively match the appar¬ 
ent magnitude distributions of RCs in Me William & Zoccali 
(2010). Thus, despite its simplicity, this model is valuable in 
understanding the kinematics of the Galactic X-shape. 

For consistency, we adopt the same solar configuration as in 
Shen et al. (2010) and Li & Shen (2012). Following G14, we 
also split the particles in a field into two samples to illustrate 
the kinematic difference between the two sides of the bar/X- 
shape. Particles within Rq from the Sun are denoted as “near”, 
and those beyond Rq are denoted as “far”. We calculated the 
heliocentric radial velocity (or line-of-sight velocity V\os) and 
proper motion (/i*, /j.^) of particles inside the bar/bulge region 
(Galactocentric radii less than 4.5 kpc). Along the minor axis 
of the Galactic bulge, some specific fields are carefully stud¬ 
ied in order to allow for detailed comparisons with observa¬ 
tions. Across the bulge region, particles are binned in the (/, b) 
plane to map the spatial trend of the mean velocity and veloc¬ 
ity dispersion. In observations, jii from multi-epoch measure¬ 
ments needs a factor cosb to represent its angular speed in 
the longitudinal direction, as required by the metric of spher¬ 
ical coordinates. In this paper we also use = jii cosb as the 
notation of the metric-corrected longitudinal proper motion. 

3. RESULTS AND DISCUSSION 
3.1. Distance Distribution towards the X-shaped Bulge 

From the solar perspective, the Galactic X-shape fea¬ 
tures double-peaked magnitude/distance distributions in fields 
passing through it (Li & Shen 2012). In this section, we study 
the distance distributions towards the X-shape in this model 
in greater detail, and investigate the properties of the double- 
peaked feature across the bulge region. 

We use kernel density estimator to test if the distance distri¬ 
bution is double-peaked (Silverman 1981; Hall & York 2001). 
In each field, the distance distribution is constructed from the 
sample particles with a Gaussian kernel. The multimodality of 
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the distance distribution is illustrated with the critical kernel 
sizes for k modes (/z^), and its significance levels are evaluated 
using bootstraping. A relatively large and statistically signifi¬ 
cant hk implies the existence of more than k modes. For each 
field, we search for the critical kernel size for a single-peaked 
or double-peaked distribution (hi and /z 2 , respectively), and 
test their significance levels (pi, P 2 ) via standard bootstrap¬ 
ping. The details of the method is described in Appendix A. 

To illustrate the properties of the double-peaked feature, 
we locate the two peaks in the distance distribution using the 
mean value of hi and /z 2 as an optimal kernel size, and find the 
separation of the two peaks in the distance distribution and the 
ratio of peak amplitudes (far to near). It is worth pointing out 
that the algorithm here is relatively conservative, while the vi¬ 
sual identification (with a priori knowledge about the double 
peaks) is more flexible in distinguishing modes. 

Figure[^shows the properties of the double-peaked distance 
distributions across the bulge region. Inside the regions where 
the double-peaked feature is robustly detected, we find some 
prominent trends. The separation of the two peaks in the dis¬ 
tance distribution increases with \b\ due to the X-shape geom¬ 
etry (upper left panel), which is especially clear in the vertical 
slice of / = 0° (upper right panel). While at the same latitude, 
the separation is roughly constant. Meanwhile, the ratio of 
peak amplitude (far to near) decreases with /. Since the bar is 
tilted, its near and far sides dominate at / > 0° and I < 0° re¬ 
spectively, as illustrated in the latitudinal slice of /? = 5° (lower 
right panel). As I decreases, the peak amplitude decreases at 
the near side but increases at the far side, which further in¬ 
creases the ratio of peak amplitude (far to near). In addition, 
at both sides, the identified peaks are closer to the observer 
as I increases (lower right panel). In other regions where the 
double-peaked feature is not robustly detected (low pi), the 
separation and peak ratio are spiky since the secondary peaks 
are likely to be spurs in the under-smoothed distribution there. 
The spatial distributions of the overdensity shown in the right 
column of Figure[2are qualitatively consistent with the recent 
OGLE-IIIRC observations in Figure 9 of Nataf et al. (2015). 

The increasing peak separation with b in this model is qual¬ 
itatively consistent with previous studies (e.g. Me William & 
Zoccali 2010; Saito et al. 2011; Li & Shen 2012; Nataf et al. 
2015). The peak separation measured in (0°,6°) (^ 900 pc 
above the disk plane at the GC) is about 1.7 kpc, which 
is similar to the result in Figure 8 of McWilliam & Zoc¬ 
cali (2010). The varying peak amplitude due to the tilted 
bar/X-shape agrees with the results in McWilliam & Zoc¬ 
cali (2010) and Saito et al. (2011), and naturally leads to 
the decreasing ratio of peak amplitudes (far to near) with I 
(lower left panel of Figure [T]). Some minor disagreements ex¬ 
ist on the latitudinal range where the double-peaked feature 
is detectable. The two split RCs in McWilliam & Zoccali 
(2010) merge around b ^ 5°, and Saito et al. (2011) found 
a similar but slightly lower limit of \b\ ^4°. The double- 
peaked density profile in Wegg & Gerhard (2013) emerges at 
400 pc above the Galactic plane, corresponding to Z? ^ 2.7° at 
Ro = 8.5 kpc. In our model, the double-peaked distance distri¬ 
bution can be detected at jZ?! ^ 3° (^ 450 pc above the disk, at 
Ro = 8.5 kpc), which is similar to the result obtained by Wegg 
& Gerhard (2013). Saito et al. (2011) obtained a longitudi¬ 
nal range of |Z| <2° for the double-peaked feature, while the 
range is latitude-dependent in our model. Towards higher lat¬ 
itudes, the longitudinal range for the double-peaked feature 
in our model becomes much wider than |Z| <2° in Saito et 
al. (2011). These minor disagreements could be due to the 


different peak-finding methods adopted. 

To further illustrate the statistical significance of the 
double-peaked feature in the distance distribution, the criti¬ 
cal kernel sizes (hi, h 2 ) and their significance levels (pi, P 2 ) 
are illustrated in Figure Closer to the disk plane, most re¬ 
gions towards the Galactic bulge feature low pi, P 2 and rela¬ 
tively small hi, Zz 2 , suggesting that the distance distributions 
are single-peaked there. In two triangular regions away from 
the disk plane (\b\ > 3°), we obtain clearly larger hi with high 
significance level, while both Zi 2 and p 2 are modest there. 
The robustly large hi but insignificant ^2 here indicates that 
in these regions, the distance distribution is strongly double- 
peaked due to the presence of the vertical X-shape. Further¬ 
more, the critical kernel size hi increases with \b\, which 
also agrees with the larger peak separation at higher latitudes. 
Since larger kernel sizes are required to over-smooth the gap 
between the two modes towards higher latitudes, the increas¬ 
ing hi with \b\ also implies the vertical X-shape. 

In practice, the identification of the X-shape can be af¬ 
fected by various uncertainties in the determination of the 
distance, such as sample contamination, photometric errors, 
non-uniform extinctions and variations in the intrinsic lumi¬ 
nosity of RCs. Moreover, the location of the peak depends 
on the peak-finding method adopted. On the other hand, this 
model may not exactly illustrate the X-shaped feature in the 
MW bulge. Despite the possible disagreements between the 
model and observations, our model is still valuable in study¬ 
ing the X-shaped MW bulge for the broad similarities with the 
observed X-shape geometry. 

3.2. Radial Velocities 

Modem spectroscopic surveys have acquired precise mea¬ 
surements of stellar radial velocities, which can provide im¬ 
portant clues to understanding the dynamics of Galactic stme- 
tures. Previous radial velocity observations have revealed the 
cylindrical rotation of the Galactic bulge, indicating the secu¬ 
lar nature of its formation history (Howard et al. 2009; Shen 
et al. 2010; Williams et al. 2011). V13 studied the radial 
velocity of RCs inside the X-shape towards (0°,-6°). They 
observed an excess of approaching stars at Vios ^ -80 km s“^ 
in the bright RCs (near side), and a similar excess of receding 
stars at Vios ^ 80 km s“^ in the faint RCs (far side). In contrast, 
inside the Plant’s window (b = -8°), the radial velocity barely 
differs between the bright and faint RCs (De Propris et al. 
2011). The numerical study of G14 related the X-shape mor¬ 
phology to a minimum in the mean difference Vios between 
the two sides. Therefore, although there have been a handful 
of kinematic studies of the bulge stars focused in exploring 
the X-shape, the kinematic imprints of the X-shape are still 
not well studied. 

First we study the radial velocity in fields along the minor 
axis of the bulge. In calculating the radial velocity, we set 
the circular velocity of the Sun to be Vb = 220 km s“\ without 
considering the solar motion with respect to the Local Stan¬ 
dard of Rest (LSR). Since the model is symmetric with re¬ 
spect to the Galactic plane, two fields with opposite latitudes 
are combined in the analysis to achieve better number statis¬ 
tics. Figurej^shows the radial velocity distributions in Z = 0°, 
b = ±2°, ±4^, ±6° and ±8°. For both sides of the bar/X- 
shape, the radial velocity spans a wider range at latitudes 
closer to the Galactic plane, which is seen in the decreasing 
velocity dispersion with \b\. In the same field, we find sys¬ 
tematically smaller velocity dispersion at the far side, indicat- 
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Figure 1. Geometric properties of the double-peaked feature in the distance distribution. The left column shows the separation (upper panel, in kpc) of the 
two peaks and their amplitude ratio (lower panel, far to near) across the bulge region. Dotted, dashed and solid contour lines represent pi = 0.5, 0.68 and 0.95, 
respectively. Color bars are scaled to highlight the trends in the region where the double-peaked feature is clear. The right column shows the peak positions in 
the I = 0° plane (upper panel) and b = 5° plane (lower panel, XY-projection). Contours display the particle density in logarithmic scale, and the black dots mark 
the identified peaks with pi > 0.68. Lines for I = ±5° and b = ±5° are also marked. In the region where the double-peaked feature is clearly identified, the peak 
separation steadily increases with \b\ (upper left), and the peak ratio (far to near) decreases with / (lower left). 


ing a decreasing velocity dispersion away from the Galactic 
disk. In all these fields, the near side has more approaching 
particles (negative Vios), while the far side has more receding 
particles (positive Vios) accordingly. The excess of approach¬ 
ing/receding particles is evident in the median velocities of 
the two sides. 

The excesses of Vios at the two sides may be subject to the 
orbital motion inside the bar structure. Figure sketches a 
naive bar-supporting orbit in this configuration. In an axisym- 
metric, unbarred disk, Vios at / ^ 0° is expected to be symmet¬ 
ric with respect to 0 km s“^ at both sides. If it has coherent 
non-circular motions, i.e. the alignment of elongated orbits 
have a preferred direction, asymmetric Vios distributions with 
respect to 0 km s“^ should be detected at both sides. In this 
case, an excess of approaching particles is expected at the near 
side, and an excess of receding particles should be observed 
at the far side, as expected from Figure Although the re¬ 
ported excesses at Vios ^ ±80 km s“^ in VI3 are not clear in 
this modej^ there are indeed more approaching or receding 
particles at the two sides, which agrees well with Figure 
Nidever et al. (2012) also used bar orbits to explain the cold 
high-velocity peaks seen in the commissioning observations 
of the Apache Point Observatory Galactic Evolution Experi- 

^ We notice that at (0°,-6°) alone, there are clear excesses at Vios ~ 
±80 km s“k which, however, become less noticeable in the combined field 
(0°, ±6°). Perhaps we need better number statistics to confirm this, both in 
observations and models. 


ment (APOGEE). However, as Li et al. (2014) demonstrated 
with numerical simulations, bars are not expected to display 
such a high velocity peak; instead, the high velocity stars are 
located around the tangential point between the bar orbits and 
the sight lines. 

Then we study the mean radial velocity (Vios), and its dis¬ 
persion (cTios) across the whole bulge region. As illustrated 
in Figure!^ Vios at both sides shows clear cylindrical rota¬ 
tion, and the longitudinal gradient of Vios is strong at jZ?! ^ 3°, 
which coincides with the latitude where the X-shape is de¬ 
tectable in this model. The zero-velocity line of Vios shifts 
towards / > 0° at the near side, and towards / < 0° at the far 
side. The shift of zero-velocity lines agrees with the coherent 
alignment of elongated orbits. As sketched in Figure]^ the ra¬ 
dial velocity is zero only close to the two tips of such orbits. If 
the bar structure is populated by such orbits, the zero-velocity 
lines at the two sides of the bar will appear at two opposite 
sides of / = 0°. Also, the decreasing radial velocity dispersion 
away from the GC indicates a kinematically hot inner bulge. 
We find that, the difference of Vios between the near and far 
sides (near - far) is generally negative, while the difference of 
(Jios is mostly positive. Since on the (/,/?) plane, the particles 
at the near side are closer to the disk plane than those at the 
far side, the positive Aaios and negative AVios may reflect the 
kinematic gradient in the vertical direction. Near \b\ = 4°, wc 
also notice the minimum of AVios along 1 = 0°. Considering 
both the absolute value and the relative depth with respect to 
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Galactic longitude [deg] Galactic longitude [deg] 

Figure 2. Testing the double-peaked distance distribution using the kernel density estimator. The left panel shows the map of hi, and right panel shows h 2 . 
Units are in kpc. Dotted, dashed and solid contour lines indicate the significance levels (pi and P 2 ) of 0.5, 0.68 and 0.95, respectively. Histograms of distances 
for several typical fields are also shown in stamps with their field positions (/, b). In two triangular regions above and below the disk plane i\b\ > 3°), we obtain 
large hi at high significance levels but plain /?2 with low p 2 , indicating that the distance distribution is strongly bimodal there. We also find that, the critical kernel 
size hi increases with \b\ near / ~ 0°, indicating that larger kernels are required to over-smooth the distribution at higher latitudes. This is consistent with the 
increasing peak separation with \b\. 






Figure 3. Radial velocity distributions at the near and far sides in fields along / = 0°, including the same field (0°,-6°) investigated in VI3. For each panel, 
two 1.5° X 1.5° fields with opposite b are combined to obtain better number statistics, and the field position {l,b) is labeled at the upper right. The near and far 
sides are shown in red solid lines and blue dashed lines respectively, with their median values (m) and dispersions (a) listed in the legend. The vertical segments 
indicate the median values of Uos at the two sides. We find more approaching (negative Vios) particles at the near side, and more receding (positive Vios) particles 
at the far side, as indicated by their negative and positive median velocities. Meanwhile, the velocity dispersions at the two sides decrease with \b\. 
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Figure 4. The face-on sketch of a bar-supporting orbit in this configuration. 
Near 1 = 0°, these orbits contribute an excess of approaching stars at the near 
side, and an excess of receding stars at the far side. Meanwhile, /j,^ at the near 
and far sides will peak at the positive and negative sides of respectively, 

where is the proper motion of the Galactic center due to the motion of 
the Sun. Its tilted alignment may also lead to the position shift in both the 
extrema of /j,^, and the zero-velocity line of Vios with respect to / = 0°. 

the GC, the minimum in our model is quite similar to G14. 
The agreement between our model and G14 indicates that the 
minimum in the mean difference of radial velocity may be 
a coherent feature of X-shaped bulges, which may originate 
from the excess of approaching/receding particles at the near 
and far sides. 

Zoccali et al. (2014) studied the kinematics of the Galactic 
bulge with the GIRAFFE Inner Bulge Survey (GIBS). With¬ 
out distinguishing the near and far sides, they mapped the line- 
of-sight velocity, and velocity dispersion of the bulge region 
using the spectra of ^ 5000 RCs observed in 24 fields. They 
found an asymmetric Vios across the Galactic longitude, and a 
vertically elongated peak of aios at the GC. Compared to the 
results in Zoccali et al. (2014), Vios of our model is nearly 
symmetric with I (left panel of Figure |^, even in the presence 
of a tilted bar. We also find a vertically elongated peak of aios 
(right panel of Figure [^. Although the velocity dispersion 
map in Zoccali et al. (2014) was interpolated from coarsely 
sampled fields, it still shares great similarity with our model. 
Therefore, the vertically elongated dispersion map in Zoccali 
et al. (2014) may be physically meaningful; it may originate 
from the velocity distribution of bar supporting orbits. We 
will study this feature in greater detail in our future work. 

3.3. Proper Motions 

The proper motion refiects the angular velocity of stars in 
the heliocentric frame. It is useful to study the dynamics 
of the MW (e.g. Kozlowski et al. 2006; Rattenbury et al. 
2007; Soto et al. 2012). Previous usage of the proper motion 
was limited by the data quality and the lack of accurate dis¬ 
tance estimation. New space missions like Gaia will obtain 
precise proper motion and distance measurements for an un¬ 
precedented number of stars, and the study of the Milky Way 


dynamics will benefit significantly from the upcoming boom 
of high-quality data. Currently, the potential of the proper 
motion is still to be fully explored, and less is known about 
the imprints of the fine structures (like the X-shape) on the 
observed proper motions. In this section, we will study the 
proper motion towards the bulge region in our model. 

First, we study the proper motion distribution along the mi¬ 
nor axis of the bulge. Figure [ 7 ] shows the proper motion in 
four fields (/ = 0°, /? = 2°, 4°,^ and 8°). To ensure robust 
statistics, as in the previous section, the fields with opposite 
latitudes are co-added to carry out the analysis. The latitudi¬ 
nal proper motion for fields at Z? < 0 is reversed according 
to the symmetry of the model. For /ij", the near and far sides 
peak at the positive and negative sides of p^iQc respectively, 

where P^iqq = Vo/Rq = -5.46 mas yr“^ is the proper motion of 
the GC due to the motion of the Sun (assuming Rq = 8.5 kpc, 
Vb = 220kms“^). The displaced distributions of /ij" at the 
two sides lead to a significant difference in their mean val¬ 
ues. For ph, the two sides have similar mean values (close to 
zero). We find a larger dispersion of and pt at the near 
side, and the dispersions at both sides decrease with \b\. The 
difference of JI^ between the near and far sides mainly re¬ 
fiects the rotation of the model bulge. In these fields, AjI^ is 
comparable to the dispersion of JI^ at both sides, which indi¬ 
cates the strong rotation of the model bulge. The role of bar¬ 
supporting orbits in this case still needs further investigation. 
In our (0°, ±6°) field we measured a relatively large AjI^ of 
2.4 mas yr“\ which is almost three times the result presented 
in VI3. Such a large difference deserves further exploration 
and explanation. 

Since the proper motion scales inversely with distance, par¬ 
ticles at the near side will show a larger proper motion dis¬ 
persion than the far side, assuming similar linear velocity dis¬ 
tributions at the two sides. This agrees well with the larger 
proper motion dispersion at the near side in this model. Here 
we test if the larger dispersion at the near side can be ex¬ 
plained solely with the shorter distance. We identify the den¬ 
sity peaks in fields along the minor axis of the bulge using the 
method described in Section 2.1, and pick the 25% particles 
closest to the identified peak position at each side. Then, we 
calculate the proper motion dispersion of these peak particles, 
and find the ratios of dispersions between the two sides (near 
to far). If the larger proper motion dispersion is mainly due 
to the shorter distance, the ratio of proper motion dispersion 
(near to far) should roughly match the ratio of peak distance 
(far to near). The results are summarized in Table 1. The ratio 
of a I is marginally consistent with, but slightly larger than, the 
ratio of peak distance, and both of them increase with latitude. 
However, o'h^Near/o'bpar poorly matches the ratio of peak dis¬ 
tance, and no clear trend with the latitude is observed. There¬ 
fore, the larger proper motion dispersion measured at the near 
side cannot be solely explained by the shorter distance. Simi¬ 
lar to the systematically smaller line-of-sight velocity disper¬ 
sion at the far side, the slightly larger (ii^Near/cripar compared 
to Rpar/RNear may be due to the decreasing linear velocity 
dispersion away from the Galactic plane. 

The mean proper motion and the corresponding dispersion 
maps across the bulge region are shown in Figure [8] Away 
from the disk plane, we find that at the same b, the extremum 
of shifts away from 1 = 0°; the maximum at the near side 
shifts towards / < 0°, while the minimum at the far side shifts 
towards / > 0°. Close to the disk plane, in the near side, the 
weak maximum of moves further away from 1 = 0°, and a 
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Figure 5. The mean line-of-sight velocity (Vios) (upper) and the velocity dispersion (cjios) (lower) across the bulge region, in km s ^ The left and middle 
columns show the near and far sides respectively, with the difference between the two sides shown in the right column (with contours). The cylindrical rotation of 
the bulge is seen in Vio^ at both sides, and cJios features a kinematically hot Galactic center. The zero-velocity line of Vw shifts away from 1 = 0°, which appears 
at / > 0° at the near side, and / < 0° at the far side. Such shift agrees with the simple bar orbits sketched in Figure]^ The minimum of AVios along I ~ 0° 
described in G14 is also clearly visible in our model (upper right). 



Figure 6. The mean radial velocity (left) and the corresponding dispersion (right) across the bulge region, in km s“^ without distinguishing the near and far 
sides of the galactic bar/bulge. The mean V\os (left), which clearly shows cylindrical rotation, is nearly symmetric with respect to / = 0°. Dashed and solid 
contours show negative and positive values, respectively. The zero-velocity line has no significant deviation away from 1 = 0°. Note that the dispersion of Vios 
features a vertically elongated peak near the GC, which agrees with the observed kinematic map in Zoccali et al. (2014). 
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Figure 7. Proper motions of the near and far sides at / = 0°, Z? = 2°, 4°, 6° and 8°. The co-added from particles at Z? < 0° are taken with the opposite sign as 
required by the symmetry of the model. Particles at the near side are in red, while the far side is shown in blue. Histograms show the distributions of and jib at 
the two sides, with the dashed lines showing their Gaussian fit, and dotted lines marking their mean values. The distributions peak at the positive side of 
(black solid lines) at the near side, and the negative side of at the far side. While for /j.^, both near and far sides peak around 0 mas yr“^. The dispersions of 
II* and fib are always larger at the near side. For both sides, all the dispersions decrease with \b\. 


Table 1 

Distance ratio and proper motion dispersion ratios between the near and far 


sides 

Latitude 

RFar / RNear 

^UNear/ ^l,Far 

^b,Near/^b,Far 

40 

1.1288 ±0.0095 

1.1426 ±0.0447 

1.1309±0.0433 

5° 

1.1881±0.0111 

1.2616 ±0.0609 

1.2253 ±0.0571 

6 ° 

1.2292 ±0.0140 

1.3629 ±0.0898 

1.1006 ±0.0788 

70 

1.2782 ±0.0155 

1.3982±0.1121 

1.3245 ±0.1248 


Notes: Rpar/RNear indicates the ratio of peak distance to the Sun (far to near). 
^i,NearFar and cTb^Near/(^b,Far indicate the dispersion ratio (near to far) for 
11*1 and lib. Uncertainties are derived using the bootstrap method. 

local minimum emerges at / > 0°. Meanwhile, at the far side, 
the minimum of /ij" shifts further away from the GC, and a 
local maximum is seen at the negative I side of the GC. 

The behavior of JI^ can be boldly explained with the toy 
model illustrated in Figure For an axisymmetric disk with 
cylindrical rotation, / = 0° is expected to be the maximum 
of at the near side, and the minimum of Jlf at the far 
side. Due to the coherent alignment of bar-supporting or¬ 
bits, the extremum of JI^ will appear close to the bar’s mi¬ 
nor axis where stellar velocity is higher. This agrees with the 
shift of the maximum/minimum regarding / ^ 0°. Close to 
the Galactic disk, the strong random motion would make the 
maximum/minimum less prominent. Compared to the higher 


latitudes where the inner buckled region of the bar dominates, 
the bar structure is longer and thinner in the disk plane. Its 
wider geometric extent leads to a larger displacement of the 
maximum/minimum close to Z? ^ 0°, and the lower velocity at 
the tip of bar orbits may result in the extrema observed at the 
other side of the maximum/minimum. 

For both near and far sides. Jib shows diagonal patterns 
across the bulge region, featuring a characteristic mean value 
of ^0.3 mas yr“^ (corresponds to ^ 12 km s“^ at the bulge 
distance). This may be caused by the projection of the cylin¬ 
drical rotation onto the inclined line of sight. For exam¬ 
ple, at the near side of (6°,6°), JI^ = -0.2masyr“\ while 
Vios = 51 km s“^ Given the mean distance of the particles 
(7.6 kpc). Jib = -0.2 mas yr“^ corresponds to a linear velocity 
of 7.2 km s“^ Their ratio, 7.2 km s^V^l km s“^ =0.141, is 
close to the tangential of the latitude (tan6° = 0.105). 

PI3 studied the proper motion of RCs inside the strip 
-1° < I < 1°,Z? ^ -5°. The near and far sides, identified 
in a probabilistic way, show an asymmetric mean difference 
about / = 0°. At the positive longitude side where the near 
side of the bar dominates, increases with I and AjI^ de¬ 
creases. While at the negative longitude side, both AjI^ and 
Ajib fiuctuate around constant values. The transition of AjI^ 
and Ajib was explained as the asymmetric streaming inside 
the X-shape. No previous model has predicted such a feature. 
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as mentioned in PI3. So we try to zoom in the same region 
in our model and check whether or not the X-shape can really 
lead to such a break of AJl^ and Ajlfy at different longitudes. 
Figure shows the mean difference of proper motion in the 
strip -1^ < I < = -5° between the near and far sides. 

Both AjI^i and are roughly symmetric with respect to 
/ = 0°; no break near I = -0.1° is seen. Therefore, it is un¬ 
clear how the bar/X-shape can really lead to the sharp break 
reported in PI3. 

We also map the proper motion dispersions (a/, at) across 


the bulge region. The first and second rows of Figure 10 


show a I and dh measured at the two sides. In these pan¬ 
els, the proper motion dispersions decrease away from the 
GC, indicating a kinematically hot central region. The far 
side in our model has systematically smaller dispersions, 
since the particles are at larger distances than the near side. 
In the field analyzed in VI3, a/ and dh for both bright 
and faint RCs are found to be around 3 mas yr“^. How¬ 
ever, our model yields smaller values in the same field, 
with (di,dh) = (2.40,2.48) mas yr“^ for the near side, and 
(1.66,1.96) mas yr“^ for the far side. This is inconsistent with 
the result in VI3 where the near and far sides share similar 
proper motion dispersions. In agreement with our model, the 
proper motion dispersions in PI3 are smaller at the far side, 
although there are large field-to-field variations. 

The distribution of proper motion in (/i*,/i/,) plane partly 
represents the velocity anisotropy in a certain sample. We 
adopt two indicators to illustrate the proper motion anisotropy 
in this model, namely the dispersion ratio, and the cross¬ 
correlation factor. The dispersion ratio (dijdi,) of a sample 
compares the random motions between the longitudinal and 
latitudinal directions. In this model, the dispersion ratio is 
saddle-shaped at both sides (the third row in Figure p^. It is 
larger than unity close to the disk plane, indicating stronger 
random motions in the longitudinal direction. In two regions 
above/below the GC, the measured ratio is smaller than 1, im¬ 
plying that the random motion may be stronger in the verti¬ 
cal direction. The proper motion distribution in VI3 is less 
anisotropic compared to our model. Both the bright and faint 
RCs in VI3 have dispersion ratios very close to 1, whereas we 
measure 0.92 and 0.84 for the near and far sides, respectively. 
The sub-fields of PI3 are spatially concentrated, but the re¬ 
sults have large field-to-field variations, making it difficult to 
directly compare our model with the results. 

The cross-correlation factor (dihjdidh, the same as the 
Pearson product-moment correlation coefficient r) measures 
the degree of anisotropy in the (/i*,///,) plane (e.g. Rattenbury 
et al. 2007). It ranges from 0 (no correlation) to ±1 (perfect 
linear correlation), where positive and negative signs repre¬ 
sent correlation and anti-correlation, respectively. However, 
it is worth pointing out that the factor here is insensitive to the 
case that the direction of anisotropy aligns with the coordi¬ 
nate axes. In this model, the cross-correlation factor features 
diagonal patterns (bottom two panels of Figure [T^. Since the 
dispersion ratio is mostly larger than 1 near the disk plane, the 
diagonal pattern implies that the direction of anisotropy tilts 
towards (but may not exactly pointing to) the GC. Towards 
higher latitudes (\b\ >3°) along / ^ 0°, the dispersion ratio 
is smaller than 1, and the cross-correlation factor is close to 
zero. Therefore, the direction of proper motion anisotropy is 
nearly vertical there. 

We also notice that at the far side, the area where a/ /a/, > 1 
is larger than the near side. Meanwhile, the cross-correlation 


factor tends to be larger at the far side. Since the particles 
at the far side have greater heights than those in the near 
side, such results indicate that the velocity anisotropy may 
be stronger at larger heights. However, the analysis here lacks 
information on the distance and the line-of-sight velocity, so 
we cannot directly illustrate the three-dimensional velocity 
anisotropy inside the bar region. We will study the velocity 
anisotropy of this model and its observational imprints in our 
future work. 

Rattenbury et al. (2007) mapped the proper motion disper¬ 
sion across the bulge region without distinguishing the two 
sides of the bar. Compared to their results, the proper motion 
dispersions are ^0.3 mas yr“^ smaller in this model (upper 
two panels in Figure [TT]). For the proper motion dispersion 
ratio, our model well matches the observed results in Ratten¬ 
bury et al. (2007). If we ignore the direction of the corre¬ 
lation, i.e. the sign of ditldidh, and compare the degree of 
anisotropy using its absolute value only, our model is slightly 
more anisotropic compared to Rattenbury et al. (2007), as 
shown in the lower right panel of Figure [Tl] 

3.4. Three-dimensional Kinematics inside the X-shape 

To explore the origin and properties of the X-shaped feature 
in the MW, one may need to study the three-dimensional ve¬ 
locity distribution inside the X-shape and then distinguish the 
possible orbital families that result in this feature. Previous 
studies suggested that the boxy/peanut bulge and the X-shape 
correspond to the banana orbits, i.e. xl orbits with the 2:1 
vertical resonance (e.g. Pfenniger & Friedli 1991; Patsis et al. 
2002; Athanassoula 2005), but other orbital families were also 
suggested to support the X-shape (e.g. Portail et al. 2015b). 
In this scenario, we would expect stronger coherent motions 
inside the X-shape due to the clustering of certain resonant or¬ 
bits. VI3 identified the streaming motion along the bar with 
the correlation of Vios vs. /i*, but they did not find any sig¬ 
nificant correlation in Vios vs. Here we study the stellar 
velocity inside the X-shape using a similar approach, and ex¬ 
amine if the X-shape possesses unique kinematic properties. 

We choose (0°, 6°) as our target field, and the opposite field 
(0°,-6°) (the VI3 field) was added (with negative pb) to in¬ 
crease number statistics. This field is chosen because the X- 
shape is well-identified there, and the sufficient density of par¬ 
ticles inside this field ensures robust statistics. For each side, 
among all the particles (“overall”), one quarter of particles 
closest to the identified peak are picked as the “peak” sample, 
while “wing” sample contains one quarter of particles furthest 
away from the GC. Rather than using a distance cut, the sam¬ 
ple selection is based on the fraction of particles away from 
the peak position, since the “thickness” of the X-shape is not 
well-defined. For each sample, we use the Pearson product- 
moment correlation coefficient (Pearson’s r) to indicate the 
significance of the correlation. We also perform linear re¬ 
gressions to find the best-fitting slopes of these samples, since 
such slopes may be weakened in the presence of uncorrelated 
data, which can then be used to compare the peak particles to 
the rest samples. 

The results are shown in Figure [T^ We find that the peak 
samples (green dots) have similar proper motion and radial 
velocity distributions with the overall samples (contours). In 
other words, the peak samples do not contain or represent any 
distinct kinematic sub-groups. However, for vs. and 
Ph vs. Vios, the slopes and Pearson’s r are larger for the peak 
sample, and smaller for the wing sample, except for pjy vs. 
at the near side, where the peak sample has the smallest slope 
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Figure 8. The mean longitudinal proper motion (top) and the mean latitudinal proper motion (bottom) across the bulge region, in mas yr“^. The left column 
shows the near side, and the middle column is the far side. The mean proper motion differences between the near and far sides are shown in the right column. 
At each latitude, the maximum/minimum of shifts away from 1 = 0°, which agrees with the coherent alignment of elongated bar orbits. Close to the Galactic 
plane, the maximum/minimum of JI^ is not prominent due to the large random motion there, and its larger displacement away from / 0° is in agreement with 

the extended thin component of the bar. The mean latitudinal proper motion shows a clear diagonal pattern at both sides. 




Figure 9. The mean differences of (left) and (right) atZ7 = -5°,-l° </<l°.In each panel, the blue lines show the mean difference, and the shaded area 
represents the uncertainties. Green dots with error bars are the measurements of Poleski et al. (2013) at similar positions. Inside the region, both AjT^ and 
have small variations, without any significant asymmetric feature like a break. In order to minimize the smoothing due to the large field-of-view, we reduce the 
field width to 0.25° here. 
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and Pearson’s r. For /i* vs. Vios, we do not find similar trends. 
The anti-correlation between and Vios seen in the model 
is in agreement with VI3, although they did not compare the 
“peak sample” to other stars in the same field. However, the 
correlation of vs. Vios, which is reported to be trivial in 
VI3, is of similar significance as ii\ vs. Vios in this model. 

The possibly stronger correlations in the peak samples indi¬ 
cate that the stellar velocity might be slightly more anisotropic 
inside the X-shape. If the X-shape is an ensemble of cer¬ 
tain resonant orbits, the slightly stronger anisotropy inside it 
seems to be a sensible expectation. We notice that the corre¬ 
lations here may match the case of downward-tilted banana 
orbits (cap-like ones from the edge-on view). The downward- 
tilted banana orbits have their abdomens and two tips above 
and below the Galactic plane, respectively. Close to their ab¬ 
domens, the positions where particles move downwards can 
explain the observed correlations. Given the clockwise rota¬ 
tion and the bar angle in this model, at the near side, a positive 
correlation between Vios and a positive correlation between 
Vios and and a negative correlation between and are 
expected; while at the far side, there should be a negative cor¬ 
relation between Vios and /i*, a negative correlation between 
Vios and jj^h and a positive correlation between and iijj. The 
results in this model qualitatively match the expected corre¬ 
lations. Moreover, vs. Vios mainly refiects the coherent 
alignment of bar orbits; it is less relevant to the vertical mo¬ 
tions. Therefore, we may not see stronger correlations of /i* 
vs. Vios in the peak samples. 

However, the explanation above may be inconclusive. The 
correlations due to the downward-tilted, if exist, may be seri¬ 
ously diluted by up ward-tilted banana orbits, and other orbits 
in the bar region. Moreover, the X-shape itself may be dom¬ 
inated by other orbital families (e.g. Portail 2015b), whose 
kinematic behavior is still poorly understood. More informa¬ 
tion is required to clarify the role of different types of orbits 
in the composition of the X-shape, such as their spatial extent, 
fraction inside the bar/bulge region, energy distribution, etc. 
It is worth pointing out that if the kinematics of the X-shape 
is dominated by the streaming motions due to banana orbits, 
different distributions are expected between the peak particles 
and the overall samples. However, the peak particles in Figure 
[T^show similar extended velocity distributions as the overall 
samples. On the other hand, if the banana orbits indeed make 
the X-shape, the density peaks we find along the line of sight 
may be in the region where particles move slower and stay 
longer, and thus diluting the velocity anisotropy. Therefore, 
although we find possible signatures of banana orbits from 
linear regressions, it is still too early to conclude that the ba¬ 
nana orbits eventually lead to the observed correlations. 


3.5. Proper Motion as a Probe of Bulge Rotation 

The longitudinal proper motions at / ^ 0° reflect the az¬ 
imuthal velocities at certain radii, and thus may be used to 
infer the bulge rotation. PI3 derived “the angular velocity 
of the Galactic bar”, i.e. the bar pattern speed with the 
proper motion measured at the two arms of the X-shape, as 
illustrated in FigureWe re-examine the PI3 method with 
our model in the following. Assuming the Galactocentric radii 
of the near and far peaks are Ri and R 2 respectively, and the 
mean azimuthal motion in the bar/bulge region corresponds 
to an angular velocity of O, the longitudinal proper motions 
at the near and far peaks are 


^ _Rxn-VQ ^ _-R2n-Vo 

p p ’ /^2,/ p _L p 

Kq-Ki Ko-\-K2 


( 1 ) 


In practice, it is easier to measure the proper motion differ¬ 
ence inside the same field. Thus, the mean difference of pf 
between the two peaks can be written 


. * _ Wl+/?2)^-(/?l+/?2)Vo 

M/ Ml,/ M2,/ (7 ?o-/?,)(7?o+/?2) 


( 2 ) 


Since R\ and R 2 are small compared to Rq, let A/? = /?i 
given pfoc = we find 


Ri 




(3) 


and AR can be approximated with the separation between the 
two peaks along the line of sight. This is essentially same as 
the equations in PI3, where they use the /-band magnitude 
difference of two red clumps (AIrc) to represent the distance 
between the two peaks (AR). 

With the procedure above, PI3 obtained a bar pattern speed 
ofQ = -87.9 km s“^ kpc~^; this value is significantly higher 
than the pattern speed of this model (38.5 km s“^ kpc“\ Mol- 
loy et al. 2015), and all current estimates of the pattern speed 
of the Galactic bar, which ranges from 30 to 60 km s“^ kpc“^ 
(Gardner & Flynn 2010; Gerhard 2011; Minchev et al. 2011; 
Wang et al. 2012, 2013; Long et al. 2013; Antoja et al. 2013). 
We derive an O of ^ 96 ± 25 km s“^ kpc“^ at (0°,5°) for our 
model, a value similar to PI3. If Q indeed represents the bar 
pattern speed, it should be nearly constant across various \b\. 

Here we further test the PI3 method by extending Equa¬ 
tion 3 to different latitudes. Table 2 shows the angular speed 
ft derived from the peak separation and Ajjf in our model. 
Here “peak particle” means one quarter of particles closest to 
the identified peak at each side. The separation between two 
peaks gradually increases towards higher latitudes, while for 
both peak particles and the overall samples, the variation of 
Ajlf is insignificant. Consequently, the derived angular speed 
ft decreases rapidly with \b\ in both samples (by nearly a fac¬ 
tor of two from |Z?| = 4° to 7°). As we have demonstrated, 
Q varies significantly with \b\, so it cannot represent the bar 
pattern speed. 

We suggest that ft (Equation 3) refiects the mean azimuthal 
rotation velocity of the bar/bulge at certain radii and heights, 
which varies with the latitude and peak separation. At higher 
latitudes, the measured Q probes the azimuthal rotation ve¬ 
locity at larger radii because the separation between the two 
peaks increases with \b\. 

To better understand this behavior, we compare the mean 
azimuthal velocity of peak particles with the azimuthally aver¬ 
aged rotational velocity of all the particles at the same Galac¬ 
tocentric radius, regardless of their height Z. We study a series 
of fields within 3.5°<|Z?|<8° (with an interval of 0.5°) along 
1 = 0°. Eor each field, the average rotation velocity (<Vs>) 
and the corresponding angular velocity (f]' =< > //?nare 

computed based on peak particles. An average radius (AR/2) 
for each pair is obtained by averaging the Galactocentric ra¬ 
dius of the two peaks. Of course, this radius increases for 
higher Galactic latitudes. Therefore, Q from Equation 3 at 
different Galactic latitudes can be naturally translated into a 
function of radius. Eigure 14 shows < > and ft' as a func¬ 

tion of the Galactocentric radius AR/2. The solid line repre¬ 
sents the azimuthally averaged < >. Apparently, < L/, > 

and ft' of peak particles (blue dots) agree very well with the 


Note that < > //? is different from the conventional angular fre¬ 

quency derived from the circular-speed curve. 
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Table 2 

Angular velocity of the Galactic bar/bulge from A/Ij" 


Lat. 

(deg) 

Separation 

(kpc) 

Afi* 

(mas yr“') 

Q 

(km kpc“^) 

AJI^ (Peak) 

(mas yr“^) 

Q (Peak) 

(km kpc“^) 

4.00 

1.07 

2.41 

121.56A37.31 

2.37 

120.14 ±37.55 

4.50 

1.28 

2.44 

107.12 ±29.96 

2.38 

105.36 ±30.05 

5.00 

1.53 

2.49 

96.38 ±24.55 

2.48 

96.07 ±24.75 

5.50 

1.63 

2.39 

89.62 ±22.71 

2.52 

92.92 ±22.90 

6.00 

1.81 

2.43 

84.65 ±20.13 

2.58 

87.89 ±20.36 

6.50 

2.02 

2.33 

76.96 ±17.81 

2.49 

80.12 ±18.04 

7.00 

2.16 

2.33 

74.09 ±16.46 

2.40 

75.36 ±16.66 


Notes: Separation indicates the distance between two peaks identified with kernel density estimator, stands for the difference of mean proper motion, and 
is the derived angular velocity. Peak is the 25% particles nearest to the identified peaks. We assume a 5% error in distance, and 10% in p,. The proper motion 
of the GC is -5.46 mas yr“^ (assuming Vq = 220 km s“^ at Rq = 8.5 kpc). 


azimuthally averaged values of all the particles (black solid 
line). Therefore, ft in Equation 3 actually measures the az¬ 
imuthal streaming motions (< 1^0 >), rather than the bar pat¬ 
tern speed. 

3.6. Locating the Galactic Center with the X-shape 

Since the X-shape is symmetric in the I ~ 0° plane, it can be 
used to locate the position of the GC (MeWilliam & Zoccali 
2010; G14). The GC can be located by fitting the peak posi¬ 
tions in the vertical I = 0° plane with the X-shape geometry. 
G14 compared two ways to locate the GC position, namely fit¬ 
ting the four arms of the X-shape, versus fitting the X-shape 
with two crossed lines. They found that Method II yields an 
accuracy better than 5% in most cases, even if the peak dis¬ 
tance has an uncertainty of 20%. However, the sample lim¬ 
itations (e.g. latitude coverage, number of fields, etc.) may 
affect the results. We try to locate the GC in our model with 
these two methods, and further beyond G14, we test their ac¬ 
curacies under different latitude coverages and total number 
of fields. 

We select N pairs of fields along the minor axis of the bulge 
(/ = 0°), with N fields above and N fields below the Galactic 
plane respectively. Fields are randomly selected in the lati¬ 
tude range where the X-shape can be clearly identified in this 
model (3° < j/?! < 8 °). In each field, the two peaks are iden¬ 
tified with the same procedure described in Section 2.1. We 
then fit the arms of the X-shape with four straight lines, and 
find their intersection points (Method I). In the ideal case, two 
of their four intersection points are right above/below the GC 
(i.e. on the vertical axis of the Galaxy), and the other two 
points are exactly on the Galactic plane. Thus the GC posi¬ 
tion can be located from the Galactic vertical axis, and the 
Galactic plane. Meanwhile, we try to locate the GC by di¬ 
rectly fitting the X-shape with two crossed lines (Method II). 
The two methods are illustrated in the left panel of Figure p7| 
We run a series of tests with different N to demonstrate the ef¬ 
fect of limited field number. To illustrate the impact of limited 
latitude coverage, in another series of tests, we also change 
the latitude range to 3° < Z? < 5.5° (“lower”) or 5.5° <b < 8 ° 
(“higher”), and compare the results with the case of full lati¬ 
tude coverage (3° <b < 8 °). 

We find that, fitting the X-shape with two crossed lines 
(Method II, right panel of Figure[T^ yields better results com¬ 
pared to fitting the four individuaTarms (Method I, middle of 
Figure p~5] ), which agrees with the conclusion in G14. Mean¬ 
while, increasing the number of fields leads to smaller uncer¬ 
tainties for both methods. We also find that, the accuracy of 
fitting depends on the latitude range of the fields used. The 


position of the GC can be better constrained with the fields 
at lower latitudes, while the fields at higher latitudes give 
the worst results. In principle, fields at higher latitudes put 
stronger constraints on the GC position. However, since the 
X-shape is well-defined at lower latitudes (3° < / < 5.5°) due 
to better number statistics (also seen in the upper right panel 
of Figure 1), the fitting results from these regions are more 
accurate compared to higher latitudes (5.5° < I < 8 °). 

Fitting the X-shaped overdensity is a novel and accurate 
way to determine the position of the Galactic center. How¬ 
ever, it still has some drawbacks. First it relies on the identifi¬ 
cation of peaks in the apparent magnitude distributions of the 
sample, which could have large uncertainties in real observa¬ 
tions due to the presence of scatters in the intrinsic luminosity, 
source incompleteness and sample contamination, etc. In the 
conversion from magnitude to distance, reddening and extinc¬ 
tion could also contribute to the uncertainties. 

Even if the uncertainties due to photometric distance de¬ 
termination can be greatly reduced with geometric distance 
measurements like Gaia does, there remains uncertainties in 
Method II. Ideally, if the peaks at the near side have similar 
heights as the far side, fitting two diagonal lines should trace 
the exact position of GC. However, due to the projection ef¬ 
fect, the peaks identified at the near side are closer to the disk 
plane than the far side. Consequently, the intersection of the 
two best-fit crossed lines tends to be closer to the Sun, lead¬ 
ing to a slight systematic underestimation of the GC distance. 
Since the projection effect is stronger at higher latitude, fields 
at lower latitudes are expected to yield better results than those 
at higher latitudes, which agrees with the right panel of Fig¬ 
ure [T^ In reality, the number of observational fields are usu¬ 
ally not symmetric above and below the Galactic plane, which 
may also affect the best-fit crossed lines to locate the GC. 

4. CONCLUSIONS 

We explore the kinematic properties of the X-shaped Milky 
Way bulge with an A^-body model that well represents the ob¬ 
served morphology and dynamics of the Galactic bulge. From 
the solar perspective, the two density peaks on the sightlines 
due to the X-shape are identified using a Gaussian kernel den¬ 
sity estimator with adaptive kernel size. In this model, the dis¬ 
tance distributions appear to be double-peaked above \b\r^ 3° 
along / = 0 °, and the increasing separation between two peaks 
with \b\ agrees with previous observations. As a result of a 
tilted bar, the height ratio of the two peaks (far to near), as 
well as their mean distance to the Sun, features a longitudinal 
gradient at higher latitudes. 

We study the line-of-sight velocity and proper motion dis- 
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Figure 10. The proper motion dispersion (tj/, a^), the dispersion ratio (tj/ / a^) and the cross-correlation factor (tj/^/across the bulge region. The first two 
rows (in mas yr“^) show the proper motion dispersion (cr/, a^) at the two sides. The third row represents the dispersion ratio (cr//f7^), and gray lines mark the 
position where = 1. The cross-correlation factor {crib/cricfb) is at the last row, with the gray contour lines for crib!criGb - 0. The left and right columns 

correspond to the near and far sides, respectively. We find a saddle-shaped dispersion ratio, and a diagonal cross-correlation factor across the bulge region, which 
indicate that close to the disk plane, the direction of proper motion anisotropy tilts towards the GC, and high above the GC, the dispersion ellipse is nearly vertical. 
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Figure 11. Proper motion dispersions across the bulge region, without distinguishing near and far sides. The upper panels represent cr/ and cr^ towards the bulge 
region (in mas yr“^), along with the dispersion ratio cr//cr/, and the cross correlation factor cr//,/cr/cr/, shown in the lower left and lower right panels respectively. 
In each panel, the results of Rattenbury et al. (2007) are marked with boxes. Compared to the same fields in Rattenbury et al. (2007), our model has slightly 
lower velocity dispersions (~ 0.3 mas yr“^ on average), but the dispersion ratio is very similar (the observed value is on average 0.3 ± 0.6% higher than the 
model value). For the cross-correlation factor, considering the degree of anisotropy only (taking its absolute value), we find on average 2.61 ± 0.58 times higher 
value in our model compared to the observed one in Rattenbury et al. (2007). 


tributions in fields {b = 2°,4°,6° and 8°) along / = 0°, with 
the opposite fields properly scaled and added to increase num¬ 
ber statistics. To illustrate the possible imprints of the bar/X- 
shape, the particles are divided into two groups (near and 
far) according to their distances to the Sun compared to R^. 
In these fields, the line-of-sight velocity distributions signif¬ 
icantly differ between the near and far sides; the near side 
shows an excess of approaching particles, while a correspond¬ 
ing excess of receding particles exists in the far side. Due 
to the rotation of the bar/bulge, the longitudinal proper mo¬ 
tions at the near and far sides peak at the positive and negative 
sides of respectively. The measured mean differences 
of between the two sides are comparable to the proper mo¬ 
tion dispersion there, indicating the fast cylindrical rotation of 
the bar/bulge. The latitudinal proper motions at the two sides 
barely differ in their mean values. We note that the differences 
of Vios and /ij" between the two sides can be qualitatively ex¬ 
plained with the coherent alignment of bar-supporting orbits. 

We test if the larger proper motion dispersion measured at 
the near side (compared to the far side) is due to the closer dis¬ 
tance. In our model, only the dispersion ratio of between 
the two sides (near to far) is marginally consistent with the 
distance ratio of the two peaks (far to near). The results here 
indicate that the systematically larger dispersion at the near 
side is not only due to the closer distance, but also contributed 
by the differences in their intrinsic velocity distributions. 

We also map the mean line-of-sight velocity, mean proper 
motion and their corresponding dispersions across the bulge 
region, both with and without distinguishing the near and far 
sides of the Galactic bar/bulge. The mean Vios at the two sides 
features strong cylindrical rotation, and the position of the 
zero-velocity line with respect to / = 0° agrees with the co¬ 
herent alignment or bar-supporting orbits. The mean differ¬ 


ence of Vios between the near and far sides shows a shallow 
minimum along / = 0°, which is suggested as the signature 
of an X-shaped bulge (G14). Without distinguishing the near 
and far sides, Vios is symmetric with respect to / = 0°, which 
barely shows the effect of the tilted bar as reported in Zoccali 
et al. (2014). For the mean /i* at the two sides, due to the 
coherent alignment of bar-supporting orbits, the extrema shift 
to the two opposite sides of / = 0° respectively. At lower lati¬ 
tudes, the displacements of the extrema are larger than that at 
higher latitudes. This may be due to the buckled inner region 
of the bar dominating at higher latitudes, while the long and 
thin component of the bar dominates at lower latitudes. In¬ 
side -1° < / < l°,Z? = -5°,the mean differences of and jj^b 
between the two sides are roughly symmetric with respect to 
1 = 0°. The breaks of AJl^ and reported in PI3 are absent 

in our model. 

For their dispersions, aios, cri and (Jb at both near and far 
sides reveal the kinematically hot central region. The disper¬ 
sion ratio is saddle-shaped at both sides, with a diagonal pat¬ 
tern of the cross-correlation factor, indicating that the proper 
motion anisotropy is tilted towards the GC. Without distin¬ 
guishing the near and far sides of the bar/X-shape, the dis¬ 
persion of Vios shows a vertically elongated peak similar to 
Zoccali et al. (2014). Meanwhile, the proper motion dis¬ 
persion is larger in the longitudinal direction, and the disper¬ 
sion ratio {cFijcFb) increases away from 1 = 0°. We also ob¬ 
serve a diagonal pattern of the cross-correlation factor. Com¬ 
pared to the proper motion measurements in Rattenbury et 
al. (2007), this model shows slightly smaller dispersion but 
stronger anisotropy. 

We also study the orbital motion inside the X-shape using 
the correlation of velocity components. In (0°,6°), for jj^b vs. 
Vios and fib vs. both the significance of the correlations 
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Figure 12. Correlation of the proper motion ii},) and radial velocity V\os at (0°, ±6°). In each panel the contours show the density of all the particles within 
the near and far sides (“overall”). Green dots mark one quarter of particles closest to the peak identified at the given side (“peak”). Linear regressions of the 
overall sample, peak particles, and one quarter of particles furthest from the GC (“wing”) from the peaks are shown in black solid lines, blue dashed lines and 
red dotted lines, respectively. The slopes and Pearson’s r-values are shown with the legend in each panel. For vs. V\os, the slope and Pearson’s r are largest in 
the peak sample. For vs. /r*, at the far side, the peak sample also has largest slope and Pearson’s r, while at the near side, the peak sample has smallest slope 
and Pearson’s r. 



Ri + R2 Separation 


Figure 13. Measuring the angular speed Q from Ajj^ and the peak separation. The gray line represents a line of sight towards the bulge region, and the yellow 
dots show the positions of the two density peaks due to the X-shape (blue shaded regions). If X-shape is symmetric about the GC, we can obtain an angular speed 
from the mean latitudinal proper motions near two peaks, and the separation between two peaks. 
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Figure 14. Mean azimuthal rotation (< >, top panel) and the correspond¬ 
ing angular velocity (^2' =< > /R, bottom panel) of the X-shaped peaks 

in the vertical / = 0° plane as a function of Gactocentric radius (~ AR/2). 
In each panel, the blue points represent the peak particles, while the black 
solid line shows the azimuthally averaged values, regardless of the particle 
height Z. For the blue points, we also labelled the corresponding latitude \b\ 
where the particles are selected. Since <V(p > decreases with increasing the 
height, at a given radius, the azimuthally (and vertically) averaged < > is 

smaller than the rotation velocity at lower |Z| but larger than that at higher |Z|. 
Due to the X-shape geometry, the Galactocentric radius and the height of the 
peak particles (the blue dots) increase with \b\, resulting in larger < > at 

smaller radius and smaller < > at larger distance compared to the average 

value (the black line). 


(Pearson’s r value) and the slope from linear regressions re¬ 
veal slightly stronger velocity anisotropy near the X-shape. 
Except for jj.t> vs. /ij" at the near side, the correlation is weak¬ 
est for the peak sample. We note that at both sides, the di¬ 
rection of the anisotropy roughly coincides with the case of 
downward-tilt banana orbits (for b > 0°). However, the un¬ 
derlying orbital family can not be clearly distinguished here. 

With this model, we extend the azimuthal rotation veloc¬ 
ity measured in PI3 with the X-shape kinematics to various 
Galactic latitudes. We find that the angular velocity decreases 
at higher Galactic latitudes, inconsistent with being the bar 
pattern speed. The angular velocity derived from the proper 
motion difference between the two sides of the X-shape actu¬ 
ally measures the azimuthal streaming motions, similar to the 
mean rotation profile < > (R). 

The X-shape is thought to be able to constrain the posi¬ 
tion of the Galactic center due to its symmetry. We test two 
methods to do this proposed in Gardner et al. (2014) with 
our model, namely fitting the four arms of the X-shape to 


find the center and diagonally joining the four arms by fitting 
the X-shape with two crossed lines. We demonstrate that a 
global-fitting yields better results than fitting the four arms in¬ 
dividually. Also, the fields at lower latitudes provide stronger 
constraints on the GC position due to better number statistics 
there. 
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Figure 15. The left panel shows an example of locating the Galactic center with the X-shape (N = 12). Circles show the identified density peaks in the I = 0° 
plane, and the red box shows the real position of the GC. One can either fit the four arms and locate the GC using the intersection points (dashed lines and green 
box, Method I), or directly fit the diagonal pattern with two crossed lines to find the GC (solid lines and blue box. Method II). The rest two panels show the results 
from these two methods. The color dots indicate the fitted positions in two directions for each N, and the error bars show the standard deviation of the offset. The 
real position is marked with horizontal dashed lines. Their colors stand for the latitude range of the fields. For each N, from left to right, green, red and blue dots 
represent the entire, higher and lower latitudes, respectively. 
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APPENDIX 


APPENDIX A: CALIBRATED SILVERMAN’S TEST 


In this section, we present a brief description of the Silverman’s test used in Section 3.1. The Silverman’s test (Silverman 1981) 
is based on the kernel density estimation of the sample. For a univariate sample Xi, ..., drawn from an unknown probability 
distribution function (PDF) /, its kernel density estimator (KDE) is 




i=l 


where v is the random variable, K is the kernel function which is assumed to be a Gaussian distribution function, and h is 
the kernel size. Here / serves as an approximation of the unknown original PDF /. It resembles a continuous version of the 
histogram, and h corresponds to the “bin size”. The kernel size h controls the amount of smoothing; larger h leads to fewer peaks 
in the KDE with the Gaussian kernel. 

To test the null hypothesis that / has k or fewer modes (Hq) against the alternative one that / has more than k modes (Hi), we 
use the critical kernel size of k modes 

hk = mf{h\f('\h) has almost^modes}. 

Large rejects the null hypothesis, indicating that there are more than k peaks. In other words, the ^-critical kernel size is the 
smallest kernel size that makes a KDE of k peaks, for a particular sample. If the sample has more than k modes, a large hk is 
required to over-smooth the KDE, and reduce the peak number to k. On the contrary, if the sample itself has less than k modes, 
only an under-smoothed KDE can result in more than k peaks, which requires a small ^-critical kernel size. 

Eor a particular ^-mode hypothesis, now we address the statistical significance of the ^-critical kernel size from the data (hk^o). 
In principle, one can draw a series of samples with equal size n from /, find their ^-critical kernel sizes (hi), and compare the 
critical kernel size from the data to those from resampling the true PDE / using the a level test: 

< hk^o 1 ^ 1 , • • • is drawn from /} > 1-a, 

where pr represents the conditional probability, and a is the level of significance. Since the peak number decreases with increasing 
the kernel size, we can address the significance level of o by accounting for the peak number of f(x;hk^o) from resampling: 

pT{hl > hk^o} = pr{/('; o) has more than k modes | 

{Xi ,... is drawn from /}. 

However, the true density / is always unknown. In order to evaluate the significance of hk^, Silverman (1981) generated the 
samples using the method in Efron (1979): 

Yi = {Xj+hkfied/^l+hlJa\ 

where Xi are sampled uniformly from the original data allowing replacement, is the variance of the sample, o is the ^-critical 
kernel size from the data, and e/ is a series of random numbers from a standard normal distribution. Here Yi actually samples a 
representative PDE /o, which is at most ^-modal, and has the same variance as the sample. Therefore, the significance level of 
hk^Q can be conservatively estimated with such a bootstraping procedure. 

Besides its conservative nature, the Silverman’s test is also not asymptotically accurate. One way to improve the accuracy of 
Silverman’s multimodality test is presented in Hall & York (2011), by using a bootstraping distribution function 

4(A) = W{k/K^ < A|4,...,Y4}. 

As suggested in Hall & York (2011), there exists a unique absolute constant Aq, such that the bootstrap test is asymptotically 
correct: 

pr{pr{%//i^,o < Xa\{Xi,...,Xn}} > l-a} 

—y pr}G(AQ;) ^ 1 — (t} = (t, 

and the constant Ac, can be calculated either analytically, or using Monte Carlo simulations. In addressing the significance level 
of the double-peaked feature (pi = l- a),we use the tabular value given in Hall & York (2011). 



